An LMI Approach to Stability Analysis of PWM 
DC-DC Buck Converters 



Alexander N. Churilov 

Department of Computer Science 
St. Petersburg Marine Technical University 
Lotsmanskaya Str. 3, 190008, St. Petersburg, Russia 
a_churilov@mail.ru 

Abstract 

The paper considers a DC-DC buck power converter employing pulse width mod- 
ulation and voltage feedback control. Global asymptotical stability of a periodical 
operating mode is examined. It is shown that the stability analysis can be reduced 
to study of a feasibility problem of semidefinite programming. 
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1 Introduction 

DC-DC buck (step down) power converters serve to convert a direct current (DC) voltage 
level to a lower DC voltage level [1]. One of the most popular types of power converters 
uses pulse width modulation (PWM). A pulse width modulator produces a train of square 
pulses with a constant frequency and with a variable duration. By controlling the duration 
one can control the output voltage. Last decade DC-DC PWM converters attracted much 
attention not only from engineers, but also from mathematicians and physicists. Their 
study lies in the mainstream of common interest to discontinuous and hybrid dynamical 
systems. 

PWM converters may be described by nonlinear functional-differential equations, their 
operating modes correspond to periodic solutions. Many publications were devoted to the 
existence of periodic or quasi-periodic solutions, to the analysis of local stability of such 
solutions, to bifurcations and ways of chaotization (see, e.g., [2-6]). However, most of 
investigations do not concern global behavior of solutions, i.e., a behavior for all possible 
initial conditions. 

There are two basic classes of mathematical models used for analysis and design of 
PWM power converters: discrete time models [7-10] and models based on averaging in 
time of the modulator's output or of the state space vector [11-15]. The applicability of 
the averaging technique suggests that the switching frequency is high when compared to 
the frequencies of continuous-time signals. 

In this paper we use another averaging approach, initiated by the work [16] and further 
developed in [17]. It is important to emphasize that our approach is not asymptotical 
or approximate, unlike the other averaging technique. However, it gives conservative 
estimates if the switching frequency is not high enough. Our method is based on ideas 
of the absolute stability theory [18] and leads to a system of linear matrix inequalities 



(LMI) [19]. This work is a natural extension of the LMI technique proposed in [20]. Some 
other methods for stability study of switched power converters via LMIs can be found 
in [21,22]. 

The main purpose of this paper is to provide a new computationally tractable proce- 
dure for stability investigations of PWM buck converters. It can be easily implemented 
with the help of a recently developed software for MATLAB and Scilab modelling systems, 
such as LMILab [23], LMITOOL [24] or SeDuMi solver [25] interfaced with YALMIP [26] 
or cvx [27]. 



2 Preliminaries 

2.1 Power Stage 

The part of a converter without a control loop is called a power stage. A conventional 
power stage of a buck converter is shown in Fig. 1 (see, e.g., [1]). Here V s is an input 
voltage, S is an on-off switch, L is an inductance, Co is a capacitance and R is a load 
resistance. Let U be a capacitor voltage (which is the output voltage at the same time) 
and %l be an inductor current. Then we get 

dU 1 . 1 TT di L 1 1 

^ = C- lL -R€^ -dT = -L U+ L Vs ' (1) 

when the switch is closed and 

dU _ 1 . 1 dh _ _l_ n ro x 

dt ~ C lL RCo ' dt ~ L U > [Z) 

when the switch is open. 

The switch is on with some constant switching period T. Let r be an on-time of the 
switch, i.e., a time when the switch is closed during the period. The value D = r/T is 
called duty ratio. 

The above described power stage consists of three parts — a voltage source, a low-pass 
filter and a load. The low-pass filter is designed to remove higher harmonics and pass 
only a constant component of the output voltage. However, in fact the output voltage 
has a small ripple. Neglecting this ripple, one obtains 

Us*^V s = DV s , 0<D<1, 

so the input voltage steps down. 

In general, the power stage of a buck converter can be described with a system of 
equations 

^ = ApXp + B p f, U = CpXp (3) 

with 

rf,\ fl, nT<t<nT + T, 

/(t) Hn -Livr n =0,1,... (4) 

0, nT + r <t < (n + 1)T, 
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Figure 1: Equivalent circuit for a basic power stage 

The signal f(t) is a train of square single-sign pulses of width r and of period T. The 
matrix A p is square, B p is a column and C p is a row. The matrix A p is Hurwitz stable, 
i.e., all its eigenvalues have negative real parts. 
E.g., for equations (1), (2) we have 



A p — 



-1/L ' 
1/Co -1/RC 



B p = 



V s /L 




C p = [0 1], xp 



Equations (3), (4) are the equations of an open loop system (without a control feed- 
back). In practice, the pulse width r is varied to control the output voltage U. Notice 
that equations (3), (4) apply to so called continuous conduction mode, which suggests 
that the inductance L is sufficiently large. 



2.2 Pulse- Width Modulation 

In DC-DC converters a single sign modulation is commonly employed, i.e. a modulator's 
output takes only two values, and 1. 

In this paper we will limit ourselves to a specific type of pulse-width modulation called 
natural or running modulation. The output of a modulator is 



/(*) 



nT <t <nT + r n , 

nT + r n <t <(n+ 1)T, 



n 



0,1,... 



(5) 



An input signal a(t) is compared with a sawtooth oscillating signal (a ramp) 

a r (t) = <7i + <7*0 - nT)/T, nT < t < (n + 1)T. 

Here a*, U\ are given positive parameters. A switching instant is determined from the 
relationship a(t) = cr r (t) (see Fig. 2). 

More accurately, the modulation law can be described as follows (see, e.g., [5]). If 
cr(nT) < (7i, then r n = 0. If 

a(nT + t) ><ti + o*r/T, < r < T, (6) 

then r n = T. In all other cases r n is the smallest value r G [0, T] to satisfy the equation 

a{nT + r) = a l + ct*t/T. (7) 
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Figure 2: A pulse-width modulation 

All the other roots of equation (7) are ignored. 

Such type of modulation is sometimes called modulation with a latch. In this case 
every switching interval contains at most one pulse, so such unpleasant effect as chattering 
cannot take place. 

Besides the above scheme, where the trailing edge of a pulse is modulated, there are 
other modulation laws with modulation of the front edge or of both edges [28]. 

2.3 Control Loop 

Here we consider a buck converter with a voltage feedback control. In most cases a closed 
loop system with a voltage mode control can be described by the equations 

dx 

— = Ax + Bf + q, a = Cx + ij, f = Ma. (8) 

(JjL 

Here A is a constant square matrix, B and q are columns, C is a row, and ip is a scalar. 
An operator M describes a pulse modulator, so functions a(t) and f(t) are the input and 
the output of the modulator, respectively. Suppose A to be Hurwitz stable, i.e., all of its 
eigenvalues lie in the open left half-plane. 

In the simplest case the control signal is defined by the formula 

a(t) = a(V ref -U(t)). (9) 

Here V re f is a constant reference signal and a is a gain. Combine (3) and (9) to obtain 
(8) with 

A = A p , B = B p , C = —aC p , ip = aV re f, q = 0. 

In a more general case two additional linear circuits are added to the control scheme, 
as shown in Fig. 3 (see, e.g., [1]). 
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Figure 3: A closed loop system 



Let state space realizations of the compensator and of the sensor gain be 



dt 



A c x c + B c £, a = C c x c + D£ 



and 



dx s 

—r- = A s x s + B S U, r] = C s x s + D S U, 
dt 



respectively, and 

Combining (10), (11) and (12) one obtains (8) with 





A p 





0" 




Bp 







A = 


B s C p 


A s 





B = 





Q = 







—B c D s C p 


—B C C S 


A c 









B c Vref_ 



(10) 

(11) 

(12) 



C=[-D C D S C P , -D C C S , C c ], ij = D c V ref . 

With the help of the change of variables x = x + A _1 q, equations (8) are transformed 
into the equations 



dx 

— = Ax + Bf, a = Cx + ij, f = Mo- 
di 



(13) 



with ip = ipo — CA l q. 

3 Operating Modes 

Operating modes of a converter correspond to periodical solutions of the closed loop 
system. As a rule, such solution is either T-periodical, or has a period multiple to T. 

The methods of finding periodical modes of a PWM system are well known. By 
straightforward calculations, one can obtain an explicit formula for the response of the 
linear part of a system to a square pulse signal. Taking into account periodicity one comes 
to a system of transcendental equations with respect to a pulse duration, which can be 
solved by numerical methods (see, e.g., [5,17]). 
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Another way to find a periodical solution is to apply the harmonic balance method [12]. 
Certainly, in this case the linear part of the system needs to have good filtering properties. 

In this work we will limit our considerations to T-periodical solutions, i.e., to solutions 
whose period coincides with the switching period. In addition, we require these solutions 
to be unsaturated, i.e., the cases f(t) = or f(t) = 1 are excluded. 

Let us formulate necessary and sufficient conditions for the existence of a T-periodical 
solution of (13). Define for any pair (r, t), < r < T, < t < T, a function 

a(r, t) = Ce At x{r) +^j + C (e At - I) A~ l B (14) 

with 

x{t) = -(1- e- AT ) 1 (I - e- Ar ) A- l B. (15) 
Here I is the identity matrix. 

Proposition 3.1 System (13) has an unsaturated T-periodical solution with a pulse du- 
ration t if and only if the following conditions are satisfied: 

< t < T; <t(t, t) = $(t); <t(t, t) > $(t) for all 0<t<r. 

Here $(t) = o x + aj/T. 
For this solution 

a(t) = a(r,t), < t < r, 

a(t) = a(r, t) + C [I - e A ^] A^B, r<t<T. 

Proposition 3.1 is proved readily by a direct computation. The next theorem provides 
computationally tractable sufficient conditions for the existence of an unsaturated T- 
periodical solution. 

Theorem 3.1 Consider system (13). Suppose that the matrix A is Hurwitz stable and 
the inequality 

(71 < 1p < (Ji + (7, + CA~ X B (16) 

is valid. Let there exist a number e > and a symmetrical matrix P such that the matrix 
inequalities 

P(A T + si) + (A + eI)P r< ~ BB\ 

PyO, CAPA T C T < 7 2 

are satisfied. 1 Here 7 = cr*/T — min{0, CB}. Then system (13) has an unsaturated 
T-periodical mode. 

The first inequality (17) is not linear in the assembly of variables P and e (it is bilinear 
in these variables). However, if e is fixed, we get an LMI in P. So, (17) can be checked by 
varying e in a loop. Obviously, the first inequality of (17) implies that the eigenvalues of 
the matrix A + el have non-positive real parts. Thus e can be varied from zero to —Re Ai, 
where Ai is an eigenvalue of A with the minimal absolute value of a real part. 

1 The character -< (^) denotes negative defmiteness (negative semi-denniteness) . The character >- (y) 
denotes positive definiteness (positive semi-dcfiniteness). The symbol T denotes transpose. 
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4 Stability Conditions 

The following statement is the main result of this paper. 

Theorem 4.1 Assume that there exists a T -periodical solution x°(t), cr°(t) of (8) and an 
estimate 

da°{t) 



dt 



;i8) 



is valid, where L\ is a positive constant. Let r° be a pulse duration for this solution. 

Let there exist a symmetrical m x m matrix H and scalars e, v such that a system of 
linear matrix inequalities 



L(H) -< R(e, u), HyO, v > 0, s>0 
is satisfied. Here L, R are square symmetrical matrices, 



(19) 



L(H) 



HA + A T H HAB HB 

B T A T H 

B T H 





R(e, v) 










3z//T 2 


n 2 CA — k 2 Ki 





(7* — TLi — v — £ 

— K 2 K 



-T\k\ 



k 2 A t C 

—K 2 Kx 
— K 2 K 

e 



T 



with k = —CB, K\ = —CAB, k 2 = T/n. 

Then for any solution x(t) of (8) the limit relationships 

lim [x(t) - x°(t)] = 0, lim r n = r° 
are satisfied. Moreover, the solution x°(t) is stable in the sense of Lyapunov. 



(20) 



Inequalities (19) can be easily verified with the help of standard software packages 
that are suitable for study an LMI feasibility problem [19]. 



5 Numerical Example 

Consider a converter with a power stage described by (1), (2), where R = 22 f2, C = 
47 nF, L = 20 mH, T = 400 /is, V s = 20 V. Let a control be described by (9) with a = 1, 
V re f = 13.5 V", <Ji = 4V. The parameter cr* will be chosen later. Then we have ip = 13.5, 
q = 0. 
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Conditions of Theorem 3.1 are satisfied for 



a* > 12.83 V. 



(21) 



Thus, under condition (21) the system has an unsaturated T-periodical mode. 

Apply MATLAB to get a response of the power stage to the pulse signal with a fixed 
duty ratio r/T. The direct modelling gives the following results: 



t/T 


0.1 


0.3 


0.5 


0.7 


0.9 


TL X 


0.15 


0.35 


0.44 


0.35 


0.15 



Thus TL\ = 0.44 can be chosen for the worst case. 

Applying Theorem 3.1 and the standard LMILab package for MATLAB gives 

<r. > 17.78 V. 



(22) 



Obviously (21) follows from (22). In result, under condition (22) the system has a globally 
stable T-periodical mode. 

If we neglect the ripple of the output voltage, the following approximate formula can 
be easily obtained: 

u°^v s aVref ~ ai . 

o~ * + aV s 

Choose cr* = 18 V, which satisfies (22). Then U° = 5 V, so the converter steps the voltage 
down from 20 V to 5 V. 



6 Proofs of the Theorems 

A Useful Lemma. We precede the proof of Theorem 3.1 with a lemma, which is a 
version of a statement proved in [29]. 

Lemma 6.1 Suppose that A is a square matrix, B is a column and e is a positive scalar. 
Then the inequality 



H(A + el) + (A T + eI)H r< HBB T H 

2e 

is equivalent to the relationship 

x T H(Ax + Bf) + e(x T Hx - 1) < 0, Vx, V/, |/| < 1. 

Since 

2x T HBf -2e- ±(x T HB) 2 f = - x T HBf - y/tej , Vx, V/, 

(24) follows from (23). 



(23) 



(24) 



Let (24) be satisfied. Obviously, (24) preserves if one replaces x for a vector ai, where 
a is an arbitrary scalar. Let us fix arbitrary x and /. Then (24) implies 

a 2 x T H(AeI) + ax T HBf - e < 0, Va. 

The latter relationship is satisfied for all a if and only if 

(x T HBf) 2 + Aex T H(A + el)x < 0. 

Since |/| < 1, we come to (23). 

Proof of Theorem 3.1: The proof generally follows the proof of Theorem 6.2 [17]. 
Let us check that under condition (16) system (13) has no saturated periodical solutions 
(equilibria). Firstly, suppose that (13) has a solution with f(t) = 0, that is r = 0. 
From the modulation law this means that er(0) < G\. Hence, from (13) one concludes 
that < G\. This contradicts the first inequality (16). Similarly, let f(t) = 1, that is 
t = T. From (14), (15) it follows that a(t) = a(T,t) = ip - CA~ X B. Then (6), (7) imply 
ip — CA~ l B > o\ + <7*, which contradicts the second inequality (16). Thus system (13) 
cannot have an equilibrium. 

Let us prove the existence of a T-periodical solution under condition (17). Define 
H = P' 1 . Then (17) imply H y 0, (23), and 

CAH- 1 A T C T < 7 2 . (25) 

Firstly, consider the case, when CB < and 7 = a*/T — CB. Let us take an ellipsoid 

E = { x : x T Hx < 1 } 

and define a Lyapunov function V(x) = x T Hx. From Lemma 6.1 it follows that 

dV 

— + 2e(V-l)<0 

along the solutions of (13). Thus if a solution starts in E, it remains in E for all subsequent 
times. Since 

m&xCAx = (CAH- 1 A T C T ) 1/2 , 
from inequality (25) it follows that 

max CAx < ajT - CB. (26) 

xeE 

Consider the operator of translation along the trajectories of (13): 

S : 5(0) h->x(T). 

As it was shown above, SE C E. Let us prove that S is continuous on E. It is easily 
seen that S is continuous provided that the pulse width tq depends continuously on x(0). 
Since r is the minimal positive root of equation (7) with n — 0, it suffices to require 

da(T - 0) (j* 
dr T ' 
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The last inequality can be rewritten as 

CAx(t ) + CBjt ajT. 

Thus continuity of S follows from (26). 

We have found that the operator S is continuous on a closed, bounded and convex 
set. Hence, it has a fixed point, which corresponds to a T-periodical solution 

Let us turn to the case, when CB > and 7 = a^/T. Consider an ellipsoid 

E = {x : (x + A~ 1 B) T H{x + A^B) < 1 } 
and a Lyapunov function V (x) — (x + A~ 1 B) T H(x + A _1 B). Obviously, 



dV 



2(x + A~ 1 B) T H[A(x + A~ l B) + B(f - 1)] 



dt 

along the solutions of (13). Arguing as above, one obtains that the ellipsoid E is invariant 
under the translation S and S is continuous on E . Hence, it has a fixed point on this 
ellipsoid. 

Proof of Theorem 4.1: The proof is based on the proofs of Theorems 3.3, 4.7 [17] 
and follows the averaging scheme proposed by A. Kh. Gelig in [16]. The idea of the 
approach is that the state vector x(t) varies much slower than the modulator's output 
f(t), so this output can be averaged in time [31]. 

Let x(t), a(t) be an arbitrary solution of (8). Consider the deviations x d (t) = x(t) — 
x°(t), a d (t) = a(t) - <j°(t), f d {t) = f(t) - f°(t), where f° = Ma°. Then we obtain 

-^ = Ax d + Bf d , a d = Cx d , (27) 



-mm 

' n 1 



(0, nT <t<nT + T T n 

sgn(r n - r°), nT + r™ n < t < nT + r^ ax , (28) 

0, nT + r™ ax <t<nT + T. 

Here r™ in = min{r n , t }, r™ ax = max{r n , r }. Recall that r° is the pulse duration for the 
given T-periodical mode. 

Define a sequence of average values 

nT+T 

1 f r -t° 



rp 

nT 

and consider a piecewise constant function 

v(t) =v n , nT <t <nT + T. 

Unlike the signal f d (t), which is modulated in width and position, the pulse signal v(t) 
is amplitude modulated. However, the average values of both signals coincide on each 
sampling interval. One has 

t 

hit) = v(t) + where u(t) = J (f d (s) - v(s)) ds. 



10 



Obviously, u{nT) — for n — 0, 1, The function u(t) may be considered as an 

averaged error of the replacement of f d (t) for v(t). In Chapter 3, [17] it was shown that 
u(t) is much less than v (t), provided that the switching period T is small. More precisely, 
\u(t)\ < T\v(t)\ for all t and 



nT+T nT+T 

J u(t) 2 dt<Y J v(t) 2 dt, n = 0,l,.... 

nT nT 

The following lemma plays the main role in the proof. 

Lemma 6.2 For any solution of (8) and for any n there exists a number t n such that 
nT <t n < nT + T and 

cr d {tn) o-* - ILi 

Inequalities (29) present a version of Lur'e sectoral constraints, which are common in 
the absolute stability theory. 

Proof of Lemma 6.2: Since the given T-periodic mode is unsaturated, the relation- 
ship 

a°(r ) =a°(T° + nT) = a 1 + a if r Q /T (30) 

is valid. Consider the following three cases. 

Case (i). Suppose that r n is defined as a minimal positive root of equation (7). Subtract 
a°(nT + r n ) = cr°(T n ) from both sides of equality (7) to obtain 

a d (nT + r n ) = a°(T°) - a (r n ) + a*v n . (31) 

Since r n = r° + Tv n , one can write 

da 

a°(r n ) - a (r°) = a°(r + Tv n ) - a (r°) = Tv n ^-(t a ), (32) 
where t s is some time instant lying between r° and r n . Then (31) and (32) imply 

V n 1 



SO 



a d (nT + T n ) da 
°* ~ ~dt 



< < 



a* + TLi ~ a d (nT + r„) ~ <r* - TL\ 

Hence condition (29) is satisfied with t n = nT + r n . 
Case (ii). Suppose that 

a{nT) < a u (33) 
so r n = and v n = —t°/T < 0. Let us estimate v n from below. In view of (30), one has 



a,v n = -o*t°/T = <7! - a°(r°). (34) 
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Then (33) and (34) imply 

<r.v n > a(nT) - a°(r°) = a d (nT) + a°(nT) - <t°(t°). 

Write 



(35) 



da 



ff u (o)- ff u (r ) = r, n ^(g, 

where t s is some time, < t s < r°. Since a°(nT) = er°(0), (35) implies 

d<7° 



> D n [a 



(t 8 ) > a d (nT). 



dt 

Then we come to (29) with t n = nT. 

Case (iii). Suppose that (6) is fulfilled, so r n = T and v n — 1 — r°/T > 0. Set r = 
and r = T in (6). Then 

a(nT) > <Ji, a(nT + T) > o\ + a*. 
Firstly, consider the case when a(nT) > <ji + d,. Then 

<7, Un = ^ + <t„ _ ^(rO) < ff ( n T) - <7°(t°) = <7 d (nT) + a°(0) - a°(r°). 
Arguing as above, we get 

< v n (a, - ^-(*«)) < cr d (nT), 

so we come to (29) with t n = nT. 

At last suppose that (6) is satisfied and a(nT) < o\ + o*. Since a(nT + T) > o\ + a*, 
there exists a number t s such that < t s < T and a(nT + r s ) = a\ + a*. It is easily seen 
that (29) is valid with t n = nT + t s . The proof of Lemma 6.2 is complete. 

Now return to the proof of Theorem 3.1. Inequality (19) can be rewritten as 



with 



L (H) 



L (H) -< R (e, v) 

HA + A T H HAB HB 

B T A T H 

B T H 



(36) 



Ro(e, v) 






-\C 





3z//T 2 

(7* 









'A T C T ' 








-«i 


e 






— K 



[CA - «i 



TXi — i/ — £ — T|k| 
Inequality (36) has less dimension than (19), however, it is not linear in e. 
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Let us make a change of variables y(t) = Xd(t) — Bu(t). Then (27) can be rewritten as 



^ = Ay + Bv + ABu, a = Cy + CBu. 
dt 



(37) 



Consider a Lyapunov function V(y) = y T Hy, where H satisfies (36). Define a quadratic 
form 



F(y,u,v) 





T 




y 




y 


u 


R (e, v) 


u 


V 




V 



Then (36) implies 



dV(y(t)) 
dt 



< F(y,u,v) - 5(\\y\\ 2 + u 2 + v 2 "" 



(38) 



along the solutions of (37). Here 5 is some small positive number. From Lemma 2 it 
follows that 

nT+T 



J F(y{t),u(t),v(t))dt<0 



nT 



for all n (see Chapter 3 [17] for details). Then (38) implies y(t) — > 0, v(t) — > as t — > +oo. 
Thus we come to (20). The Lyapunov stability can be proved as in [20]. 



7 Conclusion 

We propose an LMI based approach to finding conditions for global stability of an oper- 
ating mode of a PWM power converter. The LMI technique becomes a powerful tool for 
this type of problems when combined with well-known methods of calculation of period- 
ical solutions for PWM systems. The LMI approach leads to computationally tractable 
criteria, which can be easily implemented with the help of standard modelling software. 
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